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ABSTRACT 

We present mid-infrared Spitzer-IRS spectra of the well-known UX Orionis star VV Ser. We com- 
bine the Spitzer data with interferometric and spectroscopic data from the literature covering UV 
to submillimeter wavelengths. The full set of data are modeled by a two-dimension al axisymmetric 
Monte Carlo radiative transfer code. The model is used to test the prediction of iDullemond et "all 
( 2003) that disks around UX Orionis stars must have a self-shadowed shape, and that these disks are 
seen nearly edge-on, looking just over the edge of a puffed-up inner rim, formed roughly at the dust 
sublimation radius. We find that a single, relatively simple model is consistent with all the available 
observational constraints spanning 4 orders of magnitude in wavelength and spatial scales, providing 
strong support for this interpretation of UX Orionis stars. The mid-infrared flux as measured by 
Spitzer-IRS is declining and exhibits weak silicate emission features, consistent with a self-shadowed 
geometry. MIPS and SCUBA imaging shows that the disk has a small grain dust mass as low as 
0.8 x 1O _7 M0. The low apparent dust mass may be due to strong grain growth and settling. Further 
evidence for this is provided by the fact that the grains in the upper layers of the puffed-up inner 
rim must be small (0.01-0.4 /xm) to reproduce the colors (Rv ~ 3.6) of the extinction events, while 
the shape and strength of the mid-infrared silicate emission features indicate that grains in the outer 
disk (> 1-2 AU) are somewhat larger (0.3-3.0 /im). From the model fit, the location of the puffed-up 
inner rim is estimated to be at a dust temperature of 1500 K or at 0.7-0.8 AU for small grains. This 
is almost twice the rim radius estimated from near-infrared interferometry. Since larger (more grey) 
grains are able to penetrate closer to the star for the same dust sublimation temperature, a plausible 
interpretation of the data is that these larger grains have settled to the disk mid-plane in the puffed-up 
inner rim. A best fitting model for the inner rim in which large grains in the disk mid-plane reach 
to within 0.25 AU of the star, while small grains in the disk surface create a puffed-up inner rim at 
~ 0.7 — 0.8 AU, is able to reproduce all the data, including the near-infrared visibilities. 
Subject headings: accretion, accretion disks - circumstellar matter - stars: formation, pre-main- 
sequence - infrared: stars 



1. INTRODUCTION 

Edge-on and nearly edge-on disks around Herbig Ae 
and T Tauri stars are ideal objects for studying the struc- 
ture and composition of protoplanetary disks. In addi- 
tion to infrared emission from the disk, one has supple- 
mental information from absorption of starlight along the 
line of sight through the surface layers or the interior of 
the disk. This probing technique has been used by vari- 
ous a uthors to analyze the dust size distribution in such 
disks (IWood et al .112002 IDuchene , et, al.ll2003l lWo lf et al.l 
120031 iWatson fc Star)elfeldtll2004D and to determine the 
abundance of icy grain in the outer region s of these disks 
(|Thi et al.H2002t iPontonnidan et alJl2005T) . 

Because of the additional constraints compared with 
other disks, it is clear that edge-on disks are ideal targets 
for detailed studies involving spectroscopy and images at 
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many wavelengths. Given the many new sensitive ob- 
servational facilities, including space-based mid-infrared 
imaging and spectroscopy as well as near-infrared inter- 
ferometry, a powerful approach to analyzing these data 
is the simultaneous fitting with a multi-dimensional con- 
tinuum radiative transfer model of the disk. 

A particularly interesting sub-class of (presumably) 
nearly edge-on disks is "UX Orionis stars". For these 
objects the stellar flux on average is only marginally ex- 
tincted, but they undergo frequent, highly non-periodic 
extinction events, lasting a few days to weeks, in which 
the star dims strongly and becomes redder. A model in 
which dusty clumps in Kcplerian orbits tempora rily ob- 
scure the central star, as originally proposed bv iGrinirl 
(1988), seems the best explanation for the observational 
constraints. In particular, in the bottom of very deep 
minima, the color track eventually turns around and the 
colors become bluer. In the context of the "dust clump" 
model, this characteristic "blueing" in deep minima sig- 
nifies that the stellar photosphere has disappeared en- 
tirely from view and that scattered light from the disk 
surface begins to dominate in the optical wavebands. 

While it has been suggested that these extinction 
even ts may be due to comets passing in front of the 
star l)Oradv et al.l ll997). others have argued for a nearly 
edge-on disk scenario in which one narrowly looks over 
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the surface of the disk toward the star, where small hy- 
drodynamic perturbations temporarily pass throug h the 
line of sight l )Bertouti2000tlNatta fc Whitnevl2 000). The 
short time scales for these extinction events indicate that 
the perturbations should happe n close to the dust evap- 
oration radius. It was are: ued bv lNatta etaD ({2001^ and 
Dul lemond et all 1)2003(1 that the puffed -up inner rim of 
the d usty part of the disk (see also iDullemond et al.l 
I2001|) is likely to be the location where these extinc- 
tion events are produced. The inclination of the disk 
is such that the line of sight passes just over (or just 
through) the upper edge of the puffed-up inner rim, and 
hydrodynamic fluctuations can then vary the extinction 
of starlight along the line of sight on a time scale con- 
sistent with what is observed. Since the central star is 
typically not strongly attenuated outside a UX Orionis 
event, this scenario works only if the line of sight does 
not pass through many magnitudes of extinction in the 
outer regions of the disk. Hence, the outer parts of the 
disk cannot have a flaring g eometry for these source s. 

The conclusion reached by Dullemo nd et alJ l|2003|) was 
that all UX Orionis stars should be self-shadowed disks, 
in which the outer disk is geometrically thin enough to 
lie entirely in the shadow cast by the puffed-up inner 
rim, and therefore does not intersect a line of sight pass- 
ing just over the edge of this rim. A marginal flaring, in 
which the outer disk just barely appears above the inner- 
rim shadow, may be possible as well, as long as the result- 
ing extinction along the line of sight is marginal. It was 
shown by these authors that the shape of the spectral en- 
ergy distributions of all known UX Orionis objects indeed 
indicate that they are mostly non-flared (self-shadowed) 
disks, with some exceptions lying just on the border be- 
tween flaring and self-shadowing disks. This gives strong 
support for the "inclined disk" scenario. 

As part of the Spitzer Legacy progr am 'From Molecula r 
Cores to Protoplanetary Disks' (c2d) ijEvans et al.12 003'). 
we obtained 5.2-37.0/mi spectroscopy as well as IRAC 
and MIPS images of the well-known UX Orionis star VV 
Ser. The star was originally singled out for further study 
due to the presence of an extended (over 4') mid-infrared 
nebulosity surrounding it. This nebu losity is analyzed in 
detail in a companion paper l)Pontoppidan et aT] I2006L 
hereafter Paper II), in which it is argued that the nebu- 
losity is due to the quantum heating of Polycyclic Aro- 
matic Hydrocarbons (PAHs) as well as very small silicate 
or carbon grains. Imprinted on the nebulosity is a wedge- 
like dark band which we interpret as a shadow cast by a 
sma ll (less than a few 100 AU) nearl y edge-on disk (see 
also lPontoppidan & Dullemondlfe OOr)'). The presence of a 
disk shadow allows a relatively accurate determination of 
the inclination and position angle of the system, which is 
crucial for the interpretation of VV Ser as a UX Orionis 
star. 

In this paper we concentrate on analyzing VV Ser in 
terms of the UX Orionis phenomenon. To do this we 
create an axisymmetric radiative transfer model using 
a wide range of data to constrain the disk structure. 
This includes not only the Spitzer data, but also JCMT- 
SCUBA submillimeter imaging, near-infrared interfero- 
metric visibilities and optical light curves from the liter- 
ature. Although some aspects of an axisymmetric model 
cannot be uniquely constrained, we can make firm con- 
clusions about the disk geometry, grain sizes, inner rim 



structure, etc. More importantly, our analysis gives fur- 
ther strong evidence in favor of the scenario for UX 
Orionis stars put forward by IDullemond et~aT] <|2003|). 
In doing so, t his lends support f or th e puffed-up inner 
rim m odel of iNatta et alJ (|200 If) and Dulle mond et alJ 
l)2001[) . as well as for the interpretation of the near- 
IR to far-IR slope of the SEP in terms of flaring ver- 
sus self-shadowed disks (IDullemond fc Dominikl l2~004t 
iMeeus et al.ll200l|) . This scenario allows us to interpret 
all the observed data of VV Ser in terms of a single disk 
model with a simple geometrical interpretation. 

We construct an axisymmetric radiative transfer model 
of the entire system from 0.5 to 50,000 AU scales. The 
main objectives are to test the scenario in which the 
UX Ori phenomenon is caused by a nearly edge-on disk, 
whether a disk shadow can plausibly be projected into 
the mid-infrared PAH nebulosity (see Paper II) and fi- 
nally to use the model to constrain the structure and 
dynamics of the puffed-up inner rim. 

This method necessarily targets specific objects for de- 
tailed modeling, and one should take care when gener- 
alizing from one specific object to a class of disks. One 
central point to this study is to determine observable pre- 
dictions that can be applied to other UX Orionis stars, 
or even proto-planetary disks in general. 

In this paper we first describe the observations that 
have been collected (Sect. EJl and the characteristics of 
the source (Sect. 0). In Sect, ^the radiative transfer 
model is presented. Sect. El discusses how each observa- 
tional constraint is implemented into the model as well 
as possible sources of degeneracies. 

2. OBSERVATIONS 

The primary constraints on the model of VV Ser are 
provided by mid-infrared spectroscopy an d imaging ob- 
taine d with the Spitzer Space Telescope ((Werner et alJ 
2004). Mid-infrared spectra of VV Ser were ob- 
tained with the S pitzer Infrared Spectrometer, IRS 
lIHouck et alJ l200l using the Short-Low (SL) module 
from 5.2-14.5 /j,m, the Short-High module from 9.9- 
19.6 /im and the Long-High module from 19-37/im. The 
spectra were reduced using the pipeline version Sll.0.2 
and extracted with the c2d extraction routines. The 
SL spectra were extracted in a 12 pixel aperture. The 
background was estimated using the entire length of the 
slit and subtracted using a high-order polynomial fit to 
the measured background in order to minimize any ad- 
ditional noise from the background subtraction. The SH 
and LH spectra were extracted using full aperture ex- 
tractions. The short slits in the high resolution mod- 
ules cause any background subtraction to be highly PSF- 
model dependent. Consequently, the background was not 
subtracted from the SH and LH spectra. Since VV Ser 
is a bright source between 10 and 40 /im, we estimate 
the background to be a minor contribution, even consid- 
ering the extended emission present around the source 
(see Paper II). The orders and modules were matched 
by scaling to the shortest wavelength order. The scaling 
factors applied were 5-10%, and the absolute flux level of 
the spectrum is therefore considered accurate only to 10- 
20%. The spectrum of VV Ser has AOR Key 0005651200. 
We present IRAC and MIPS imaging at 3.6, 4.5, 5.6, 8.0, 
24 and 70 /zm of the area surrounding VV Ser in Paper 
II. 
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TABLE 1 
Photometry of VV Ser 



Wavelength (/im) Flux [Jy] 1 

"02 (4 ±2) x 1CT 4 

0.365 0.011 ±0.0005 

0.44 0.035 ± 0.002 

0.55 0.077 ±0.004 

0.70 0.14 ±0.007 

0.90 0.28 ±0.014 

1.235 0.47 ±0.02 

1.662 0.97 ±0.05 

2.159 1.87 ±0.09 

3.79 3.0 ±0.1 

4.64 3.4 ±0.2 

5.5 4.2 ±0.1 

8.0 4.5 ±0.1 

12.0 4.2 ±0.2 

25.0 3.4 ±0.2 

36.0 2.7 ±0.2 

70.0 0.35 ±0.1 

450 < 0.18(3o-) 

850 < 0.03(3o-) 



Reference 



IUE 

(Rostopchina & Grinin 2001) 
(Rostopchina & Grinin 2001) 
(Rostopchina & Grinin 2001) 
(Rostopchina fc Grinin 2001'! 
(Rostopchina fc Grinin 20011 
2MASS 
2MASS 
2MASS 

(Berrilli et al. 19921 
fBerrilli et al. 19921 
Spitzer-IRS, this paper 
Spitzer-IRS, this paper 
Spitzer-IRS, this paper 
Spitzer-IRS, this paper 
Spitzer-IRS, this paper 
Spitzer-MIPS, paper II 
JCMT-SCUBA, this paper 
JCMT-SCUBA, this paper 



The optical photometry is for the "quiescent" state. 

In support of the Spitzer observations, additional con- 
tinuum data at 850 and 450 /im were obtained in service 
mode with the Submillimetre Common User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope 
(JCMT) 6 on 2005 January 27. VV Ser was observed in 
a 64 point jiggle map (2.3' field) with approximately one 
hour of integration. The weather was good with 225 Gz 
sky opacities of about 0.05. The pointing was checked 
regularly and found to be accurate to within a few arc- 
seconds. Maps of CRL2688 and Mars were used for cal- 
ibration purposes. The absolute calibration is accurate 
to about ±20%. The beam size (HPBW) of the SCUBA 
observations is approximately 14" at 850 /im. No con- 
tinuum emission was detected toward VV Ser with a 3ct 
upper limit of 10 mJy beam -1 . 

An ultraviolet spectrum obtained with the Inter- 
national Ultraviolet Explorer (IUE) was taken from 
the Multimission Archive at STScI (MAST) 7 . More 
than 12 years of photometric monitoring in the op- 
tical UBVRI bands was tak en from the catalogue of 
iRostopchina fc Grininl l)200lH . The photometry for the 
VV Ser point source is summarized in Tabled 

3. SOURCE DESCRIPTION 

VV Ser is a typical UX Orionis star located in the 
Serpens molecular cloud. The adopted distance of VV 
Ser is ~ 26 pc as determined for the Serpens molecular 
complex bv iStraizvs eta!] (dH). There is some varia- 
tion in the literature regarding the effective temperature 
of VV Ser, although most classify the star as a late B- 

6 The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the parent organizations: the Particle Physics 
and Astronomy Research Council in the United Kingdom, the Na- 
tional Research Council of Canada and the Netherlands Organiza- 
tion for Scientific Research. 

7 Some of the data presented in this paper were obtained from 
the Multimission Archive at the Space Telescope Science Institute 
(MAST). STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
Support for MAST for non-HST data is provided by the NASA 
Office of Space Science via grant NAG5-7584 and by other grants 
and contracts. 



type. iHernandez et all <|2004[) identify VV Ser as a B6 
star due to the presenc e of a number of H e lines in their 
low-resolution spectra. iMora et alJ J2001) find a spectral 
type of AO using high resolution spectroscopy. The exact 
choice does have measurable consequences for the model 
SED at wavelengths below ~ 1 /im. The optical colors 
favor a star closer to AO than B6. The stellar luminosity 
of 49 Lq also suggests that VV Ser is of later type since 
a B6 cla ssification puts it below the ZAMS in the HR 
diagram <)Siess et al]l200C|) . 

Since the optical-UV spectrum of VV Ser does not 
appear to be strongly veiled and is of high quality, it 
is possible to derive an accurate stellar luminosity of 
(4 9 ±5) x (d/260pc) 2 Lq . Using the evolutionary tracks 
of iSiess et all (|2000f) and an effective temperature of 
10,200K, this gives a stellar mass of 2.6 ± 0.2 Mq and 
an age of 3.5 ± 0.5 Myr for a solar metallicity. The un- 
certainties reflect the uncertainties in the luminosity and 
T e ff rather than those of the model tracks. 

The optical colors of the star correspond to a steady 
extinction of Ay ~ 3, but with frequent, non-periodical 
dips lasting of order 10 days with brightness minima cor- 
responding to 0.5-4.0 mag of additional extinction. 

4. MODEL 

To model the observed SED and Spitzer imaging of 
VV Ser, we use the axis ymmetric Monte Carlo radiat ive 
transfer code RADMC (jDullemond k. Dominil3l2004|) in 
combination with the raytracer of the mo re general code 
RADICAL l|Dullemond fc Turollal l20"of1) . The density 
structure is axisymmetric, but the photons are followed 
in all three dimensions. The code is used to derive the 
temperature structure of a given dust distribution. The 
dust temperature is determined for a passive disk, i.e. 
it is assumed that accretion heating is negligible rela- 
tive to direct stellar irradiation. This is justified by 
the high luminosity of the central star (49 Lq), since 
the accretion luminosity is < 3 Lq for accretion rates 
M < 10- 7 M o year- 1 i|Kenvon et aUll993j) . The low 
mass of the disk [8 x 1O~ 6 M x (-R gas todust/100)], w here 
-Rgastodust is the gas to dust mass ratio (see. Sect. I5.1fl 
also limits the possible range of accretion rates, since 
only a small fraction of the entire disk mass can be ex- 
pected to accrete per year. The low mass of the disk is 
particularly well constrained by the 70 /mi MIPS point 
as well as the upper limit to the 850 /im flux. 

Once a temperature structure has been determined, 
RADICAL can be used to create SEDs and images. The 
code is restricted to isotropic scattering. Taking non- 
isotropic scattering into account may change the opti- 
cal colors of the model somewhat, but without resolved 
optical images of the disk, non-isotropic scattering ef- 
fects cannot be properly constrained. One advantage of 
the code is the ability to include an arbitrary number 
of dust components, each with a unique opacity. The 
setup has been used to mo d el similar protostellar disks 
llPontoppidan et alJ 120051: iPontoppidan fc Dullemondl 
2005), while the code itself has been extensively tested 
for other applications l^ an Bemmel DullemondJ 20031: 



iDullemond fc Dominikl 120051: iMeiierink et alJl2005[) and 

has recently been benchmarked relati ve to indepen- 
dent continuum radiative transfer codes (|Pascucci et alJ 
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4.1. Disk geometry 

VV Ser is modeled using a density structure consisting 
of a central disk surrounded by a spherically symmet- 
ric envelope. For the central disk, the following density 
structure is adopted: 



p(R,Z) = — §pL=,. N p 



H p (R)V2n 



2H v {Rf 



(1) 



where = Sdisk x (-R/-Rdisk) p is the surface density 

and 



H P (R)/R = (H disk /R disk ) x (R/Rdisk) 



(2) 



is the disk scale height with an the flaring index. The 
flaring index is 2/ 7 for a passive irradiated disk with a 
grey dust opacity l|Chiang fc Goldreichl ll997L For aa < 
the disk is non- flaring and shadowed. We have chosen 
an intermediate value of 1/7. However, once the disk is 
self-shadowed, the exact value of the flaring parameter 
does not have a strong influence on the SED. Eq. ^is 
appropriate for the vertical structure of an isothermal 
disk in hydrostatic equilibrium and gives a convinicnt 
parametrization, even for disks that are non isothermal 
and not in hydrostatic equilibrium. 

Indeed, we do not require a vertical structure deter- 
mi ned by hydrostatic eq ui libriu m such as that described 
in IChiang fc GoldreichI l|1997[) and IDullemond et all 
(|2001^l . but let the outer vertical scale height be 
a free parameter. Wh ile the unique solution of 
IChiang fc GoldreichI 1(19971) assumes that gas and dust 
are well-mixed and thermally coupled, there are several 
mechanisms that may operate to cause the vertical struc- 
ture of the dust to deviate significantly from hydrostatic 
equilibrium. For instance, dust grains may decouple from 
the gas by grain growth and settling in which case hy- 
drostatic equilibrium no longer applies. Recent modeling 
of dust settling indicates that it is a rapid process in cir- 
cumstellar disks and produce s highly observable effects 
(Dull emond fc Dominikll2004|) . 

For the structure of the inner disk, we adopt a puffed- 
up inn er rim model similar to that of IDullemond et ahl 
(2001). The scale height of the inner rim, H Ilm , is in- 
cluded as an adjustable parameter. The radius, i? r im, of 
the inner rim is set to the location at which a temper- 
ature of 1500 K is reached. Note that since the dust in 
the inner rim does not have a grey opacity, the radius 
of the inner rim is pushed to a significantly larger radius 
than that of grey dust. The puffed-up inner rim is con- 
nected to the "normal" disk structure by changing the 
scale height power law of Eq. J3J) at a few times the in- 
ner rim radius such that the scale height rises inwards to 
meet the puffed-up inner rim. The radius of the break 
in the power law is not strongly constrained, and a pre- 
sumably reasonable value of 2.5 x R- ln is adopted here. 

The disk is surrounded by a large "envelope" simulat- 
ing the presence of the large-scale molecular cloud. The 
column density of the cloud is adjusted to reproduce the 
measured extinction toward W Ser outside of extinc- 
tion events. A cavity centered on the star+disk system 
is carved out of the cloud. Paper II discusses the model 
fit on large spatial scales, including the detailed structure 
of the envelope, and how a local model grid for the en- 
velope parameters was calculated, while this paper deals 



with the parameters of the disk itself. 

4.2. Dust opacities 

The mid-infrared spectrum and the optical colors in- 
dicate that several different dust opacities must be in- 
cluded in the model, corresponding to different parts 
of the disk. The silicate emission band at 9.7 fiva. is 
broader and flatter than an inters tellar feature (see Sect. 
15. 51 1. indicative of grain g rowth l|Bouwman et alJ 120011: 
iKessler-Silacci et all200fi|) . whereas the colors of the flux 
variations at optical wavelengths indicate that the grains 
in the innermost parts of the disk are not very different 
from small interstellar grains. Grains large enough to 
produce a flat silicate feature will have grey opacities 
in the optical wavebands. Finally, the PAHs necessary 
to model the large nebulosity surrounding VV Ser con- 
tribute an opacity of their own (see paper II for details). 

The model allows for an arbitrary number of distinct 
dust reservoirs with relative abundances that may vary 
throughout the grid. In the VV Ser disk, we use two dis- 
tinct dust components: The first consists of small grains 
that produce non-grey opacities in the optical wave- 
bands, the second component consists of larger ~ 1/im 
grains. The small-grain dust is located in the puffed-up 
rim (< 1 AU), while the rest of the disk is populated with 
the larger grain component. 

The opacities of the thermalized dust grains are calcu- 
lated using Mie theory. The grains are spherical silicates 
with inclusions of carbonaceous material. The silicate 
optical constants are those o f oxygen-rich silicates from 
lOssenkopf fc Hennind l(1994lL while the carbon optical 
constants a r e of c arbon clusters formed at 800 K from 
Pager et alJ l|T998). We do not include any ice compo- 
nent in the dust opacity. The Maxwell-Garnett e ffec- 
tive medium formula (e.g. IBohren fc Huffmanll 1983(1 has 
been used to calculate the optical constants of the sil- 
icate with carbon inclusions. The volume fraction of 
the carbo n is 30%. roug hly consistent with that deter- 
mined by Drainc (2003). The final opacities are then 
calculated with a grain size distribution consisting of a 
power law: dN(a)/da oc a~ 3 5 with minimal (a min ) and 
maximal (a max ) grain sizes. The disk SED is attenuated 
by foreground dust from an extended envelope, which 
is discussed in detail in Paper II. The exact choices of 
grain size distributions for the disk dust are justified in 
Sects. 15.41 (optical colors) and 15.51 (silicate emission fea- 
ture), and the dust parameters are summarized in Ta- 
ble [21 For comparison to other dust opacities, it is use- 
ful to note that the "large grain" dust opacity of the 
disk has Kmoum = 0.01 cm 2 g _1 x (i? gast odust/100), where 
-Rgastodust is the gas to dust ratio. This is similar to, but 
on the small side of the o paciti es (for coagulated grains) 
of lOssenkopf fc Hennind" 111991 . 

5. OBSERVED FEATURES EXPLAINED BY THE MODEL 

5.1. Constraining the model 

Due to the long computing times for a single model, 
the model parameters are varied by hand. Rather than 
calculating a comprehensive grid in all parameters, hun- 
dreds of models were calculated until a good fit was found 
to all the observable quantities. While this is not an 
optimal method for quantifying any degeneracies in the 
fit, some parameters may be independently constrained 
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on physical grounds. In Table the parameters of the 
best-fitting model arc summarized. This model has been 
constructed using all available observational constraints 
apart from the near-infrared interferometry visibilities 
discussed in Sect. 15.61 We will refer to this model as 
"Model 1". The density and temperature structure of 
Model 1 are shown in Fig. □ The SED of Model 1 is 
compared to the observed SED of VV Ser in Fig. [21 In 
Sect. 15.2115.51 we discuss how each observation property 
is connected to the model. In Sect. 15.61 it is discussed 
how the near-infrared interferometry from the literature 
requires Model 1 to be modified on small scales. 

5.2. Extinction events 

The optical light curve of VV Ser, observed by 
iRostopchina fc Grininl l)2001[) . is shown in Fig. It 
is seen that the light curve of VV Ser contains a number 
of extinction events, each lasting a few days. Although 
each event is undersampled, their average duration can 
be estimated by searching for the duration most compat- 
ible with all the observed events. Assuming that the ex- 
tinction events have an identical Gaussian shape in the 
light curve and vary only in intensity, it is found that 
the events are consistently shorter than 5.5 days and 
longer than 3.0 days (FWHM). The best-fitting Gaus- 
sian is shown in the right panel of Fig. [3J The width of 
any perturbation of the inner rim causing an extinction 
event is likely to be less than or comparable to the pres- 
sure scale height of the disk. Scenarios where this is not 
the case are interesting to consider, but probably require 
hydrodynamical simulations to constrain. Assuming Ke- 
plerian rotation, this requirement constrains the radius 
of the inner rim of the disk, since a rim radius that is 
too small will cause extinction events to be shorter than 
those observed. An inner radius corresponding to a dust 
sublimation temperature of 1500 K for small grains is lo- 
cated at 0.7-0.8 AU at the VV Ser luminosity of ~ 49 L Q . 
This corresponds to a rotation period of the inner rim 
of ~ 180 days and a physical size of the dust pertur- 
bation causing the extinction in the azimuthal direction 
of 0.09-0.15 AU. A spherical dust perturbation then has 
an average density of 2-3 x 10 9 cm -3 in order to produce 
a maximum observed extinction of 5 magnitudes in the 
V band. These values are close to those of the model 
with an inner rim scale height of H X i m / R rim = 0.105 
and a maximum gas density along the line of sight of 
~ 10 9 cm -3 between extinction events for an inclination 
of ~ 70° . A specific extinction event then corresponds to 
creating an over-density in the inner rim enhanced by a 
factor of a few compared to the "quiescent" state of the 
inner disk. 

5.3. The outer radius 

The outer radius of the disk is not well constrained 
and is somewhat degenerate with the disk opening an- 
gle. The mass is constrained by the 70 /zm MIPS point 
and the 850 /xm upper limit to be as low as 8 x 10~ 6 x 
(-Rgastodust/100) Mq. As usual, it must be stressed that 
the observations do not probe grains that have grown to 
sizes of more than a few mm, and the disk may contain 
a significant population of such grains in the mid-plane. 
Similarly, any total mass estimate will necessarily assume 
a dust-to-gas ratio of 100, but this is likely to be very 
different for a disk as evolved as VV Ser. The mass of 




R [AU] 

Fig. 1. — Physical structure of the disk model. Top panel: The 
temperature of the disk as calculated by RADMC in Kelvin. Bot- 
tom panel: The density structure of the disk in units of hydrogen 
molecules per cm 3 . The three lines indicate the path from the cen- 
tral star to an observer for different inclination angles. Note that 
for the best-fitting inclination angle of 71.5°, the line of sight to the 
central star passes through both the puffed-up inner rim as well as 
part of the outer disk. 



TABLE 2 
Adopted dust components 



Component a m i n [/jm] 


Imax 


[firn] Disk regime 


Model 1 (Standard model) and Model 2 


Small grains 0.005 


0.4 


R < 1 AU, R > 5000 AU 


Large grains 0.3 


3.0 


5000 > R > 1 AU 


Model 3 


Small grains 0.005 


0.4 


0.7 < R < 1 AU, R > 5000 AU 


Large grains 0.3 


3.0 


R < 0.7, 5000 > R > 1 AU 



small dust grains is accurate to the uncertainty in the 
far-infrared opacity, i.e. by a factor of a few. 



5.4. Optical photometry and variability 

As discussed in Sect. 14.21 the observation that the 
grains responsible for the silicate emission features ap- 
pear to be different from those responsible for the extinc- 
tion events is taken into account in the model. While 
many UX Orionis stars appear to have disks in which 
the dust causing the extinctio n events is dominate d 
by small grains (a < 1/jm) llvan den Anckerl fl9991. 
models of grain growth predict that the inner parts of 
disks very quickly become d ominated by large grains 
ijDullemond fc D ominik1 l2005|) . A priori, it is not clear 
why the grains responsible for the extinction events have 
a different size distribution than those of the outer disk. 
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TABLE 3 
Best-fitting model parameters 



Parameter Model 1 



Model 2 Model 3 



Estimated range (see text) 



Lum. 
Sp. T. 
Dust mass 



49 M Q 

B9 (10 200 K) 
0.8 x 1O- 7 M 
-1 
1/7 



H disk /R disk 0.125 

tfrim/Rrim 0.105 

T sub 1500 K 2000 

R lim 0.7 AU 0.4 

i?disk 50 AU 

Incl. 71.5° 40° 

P.A. 11° 170° 



45-55 
A2-B6 

0.1-8.0xlO- 7 M Q 
not explored 
not explored 
0.1-0.15 
0.1-0.15 

1500 1300-1600 K 

0.27 and 0.80 given by T sub 
20-200 AU 

68° 65-75° 

15° 10-20° (see Paper II) 




1.0 10.0 100.0 1000.0 

Wavelength [micron] 

Fig. 2.— The observed SED of VV Ser compared to SEDs of 
Model 1. The three model SED curves correspond to a varying 
optical depth of the perturbation of the inner rim responsible for 
the observed UX Ori extinction events (the solid curve corresponds 
to AAy = 0, the dotted curve to AAy = 0.75 and the dashed 
curve to AAy = 4.0 mag.). The sets of optical data corresponding 
to the symbols (+), (x) and (□) are single points in the light curve 
obtained at JD2450624.358, JD2450683.290 and JD2450956.499, 
respectively. 

As a first attempt, we simply model the observations by 
changing the opacity of the dust in the puffed-up inner 
rim at radii between the rim and the break in the scale 
height power law defined in Sect. 0] (see Table [5J. This 
has the effect of converting the grey (in the UV-NIR) 
disk opacity to one reproducing the optical colors. Note 
that this also causes the inner rim to absorb visible pho- 
tons more efficiently than it re-emits the energy in the 
near-infrared. This significantly increases the tempera- 
ture of the star-facing surface of the rim and has the effect 
of pushing the rim outwards, given a certain dust subli- 
mation temperature. Possible alternative geometries are 
discussed in Sect. 15.61 

The grain model of the material causing the extinction 
events is constrained by the optical extinction curve that 
can be derived by the observed optical reddening vectors. 
The extinction events correspond to grains that appear 
to have properties similar to those of interstellar grains. 
The reddening vectors in the optical color-magnitude di- 
agrams are fitted by a power-law grain size distribution 
by varying the minimum and maximum sizes as well as 
the power law index. Not unexpectedly, a good fit is 
found for a m in = 0.05 /im, a max = 0.4 /im, giving an op- 
tical extinction law with Ry ~ 3.6. 

As discussed in Sect. 15. li the perturbations of the inner 
rim are presumably not wider in the azimuthal direction 
than the pressure scale height. In essence this means 



that a dust perturbation in a Keplerian orbit roughly 
subtending 10° as seen from the star is responsible for 
an extinction event. The optical variability of VV Ser is 
therefore fundamentally a three-dimensional effect and to 
model it within the axisymmetric framework presented 
here requires some assumptions. If it is assumed that 
the azimuthal angle the dust perturbation subtends as 
seen from the star is <C 2-7T, and that the shadowing and 
scattering effects caused by the perturbation will not af- 
fect the thermal part of the SED. In this case, only the 
extinction at a specific inclination through the perturba- 
tion needs to be taken into account. This means that the 
temperature structure of the system can be calculated 
for the "quiescent" state of the inner rim in which no 
localized perturbation of dust is providing extra extinc- 
tion toward the star. The "active" state corresponding 
to an extinction event is simulated by using the temper- 
ature structure calculated by the Monte Carlo code, but 
adding extra extinction to the input spectrum of the cen- 
tral star when calculating the final SED by raytracing. 
Thus, by attenuating the star by a variable amount of 
dust, model tracks through the optical color-magnitude 
diagrams can be constructed and compared to the ob- 
served tracks. The effect on the SED of attenuating the 
input stellar spectrum by a variable extinction (AAy) 
is shown in Fig. [21 and compared to photometric data 
obtained during and outside of extinction events. 

In the model for VV Ser, a blueing effect occurs nat- 
urally because scattering dominates the optical flux for 
high optical depths through the extincting blob of dust. 
The resulting model tracks are plotted on the observed 
photometric points for VV Ser in Fig. 0] The tracks cor- 
respond to Ay's of 0.0-4.0 mag. It is seen that a strong 
blueing effect appears for Ay > 2.5 mag. The blueing is 
strongest in the B — V color. In R— I a slight reddening 
is seen instead. This is due to the / band photons being 
slightly dominated by thermal emission from the 1500 
K inner rim rather than from scattering of photosphcric 
photons. The model approximately reproduces this be- 
haviour, although some differences are noticeable. The 
use of isotropic scattering, as opposed to non-isotropic 
scattering that tends to be forward-throwing, may ac- 
count for some (up to a factor of 2 in absolute flux) of this 
difference. There are also slight absolute offsets between 
the model colors and the observed colors. This offset may 
be due to inaccuracies in the stellar spectrum. The dif- 
ferences in colors between the possible range of spectral 
types (A2 to B6, see Table|3J) is 0.55, 0.20, 0.14 and, 0.12 
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Fig. 4. — Color-magnitude diagrams of the photometry of VV Ser from Rostopchina & Grinin (2001) compared to curves calculated 
using the radiative transfer model. The model curves correspond to extinction events of Ay = — 4 mag. The small 0.15 mag offset in the 
R-I color is due to an inaccuracy in either to model photospheric color or the inner rim color. 



magnitudes for U — B, B — V, V — R and R — I, respec- 
tively, differences similar to or larger than the model- 
data offsets. Finally, small changes in the Johnson filter 
curves used in the models can produce similar offsets. 

5.5. Mid-infrared spectrum 

The 5.2-37.0/mi spectrum as observed by Spitzcr-IRS 
is dominated by the silicate emission bands at 9.7 and 
18 /im. The 9.7 /urn band is significantly broader and 
flatter than that of interstellar silicate grains. This is 
illustrated in Fig. 03 where the shapes of the observed 
silicate bands are compared to those of dust opacities 
corresponding to grains of different sizes. In the fig- 
ure, the disk continuum has been subracted by fitting a 
power law to the spectr um at 5.5, 13 and 30-40 ^m, fol- 
lowing the procedure of iKessler-Silacci et alJ l)2006f) . A 
broadening of the 9.7 /im silicate band is indicative of a 
grain size distribution dominated by grains larger than ~ 
1/xm Ijvan Boekel et al.ll2005l IKessler-Silacci et al.l f2006). 
Other emission features may appear to broaden the sili- 
cate bands. These include emission bands due to PAHs 
that in the case of VV Ser are visible at 6.2 /im and possi- 
bly at 11.3/zm, although this feature may also be due to 



crystalline silicates (see also Geers et al., accepted). Be- 
cause they are bright and unresolved, the PAH emission 
features are not related to the larger nebulosity discussed 
in Paper II and are most likely due to material associ- 
ated with the disk. The exact dust size distribution of 
the "large grain" component is not strongly constrained 
by the silicate features, in particular not in the disk mid- 
plane. Here, we adopt a dust mixture with diameters 
flmin = 0.3 /im and a max = 3.0 /im. The most important 
property of the large grain mixture, for the purposes of 
this model, is that it has grey optical and near-infrared 
opacities. As long as this condition is satisfied, the con- 
clusions reached from the modeling are not sensitive to 
the details in the size distribution of the large grain com- 
ponent. It is likely that the dust mass of the disk has a 
significant population of much larger (m illimeter-sized) 
grains (e.g. Dullcmond & Dominik 2005), but this com- 
ponent is not constrained by the observations presented 
here. 

Another parameter constrained by the Spitzer-IRS 
spectrum, central to the discussion of VV Ser as a UX 
Orionis class Herbig Ae star, is the slope of the mid - 
infrared SED. As suggested bv lDullemond et all |2003), 
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Fig. 5. — The Spitzer spectrum of the silicate emission features 
of VV Ser compared to opacities of spherical grains with diameters 
of 1.0 and 3.0 fim. A power-law continuum has been subtracted 
from both the observed spectrum as well the opacity. 

if the inner-most regions of the disk are responsible for 
the extinction events, the outer disk must have a scale 
height, H/R, similar to or smaller than that of the 
puffed-up inner rim. This is the definition of a self- 
shadowed disk, characterized by a declining spectrum for 
A > 5 /im with weak silicate emission features. This type 
of object also c orresponds to a "gr oup II" object in the 
classification of iMeeus et alJ l)2001|) . VV Ser clearly has 
a mid-infrared SED that suggests that the disk is self- 
shadowed. However, in the model, the declining SED is 
in part caused by the low mass of the disk, i.e. to avoid 
predicting too much flux at A > 5 /im, the models are 
constrained both to be self-shadowed and to have a to- 
tal gas + dust mass of M < 10~ 5 Afo. In particular, 
the low mass is necessary to keep the SED declining be- 
yond ~ 30 /im, while the shadow of the puffed-inner rim 
ensures that the SED declines between 5 and 30 /im. A 
simple test of this was conducted by removing the puffed- 
up inner rim from the model, which resulted in a flat SED 
below 30 /im. In the best-fitting model, 90% (1.3 mag) 
of the extinction at 0.2 /im in the disk along a sight-line 
inclined at 71.5°relative to the disk axis occurs within 
1.5 AU. 

Does the model suggest that the scale height of the 
disk has been lowered due to dust settling, or has only 
the apparent dust mass been lowered? The scale height, 
Hp/R for hydrostatic equilibrium at 50 AU is 0.09. This 
is somewhat smaller than the model scale height of 0.12. 
However, this value is not strongly constrained, as long 
as the inner disk shadows the outer, it is difficult to tell 
the difference without having more sensitive photometry 
at wavelengths longer than 70 /im. 

5.6. Near-infrared interferometry 

Several interferometric measurements of the near- 
infrared morphology of the innermost par t of the VV 
Ser system are available in the literature. lEisner et al.l 
l)2003l l2004|) presented -fT-band interferometry observa- 
tions of VV Ser performed with the Palomar testbed in- 
terferometer. These authors find a K -band size of VV 
Ser of ~ 3 mas fo r a generic flared disk w ith a puffed- 
up inner rim from lDullemond et alJ (|2001|) . This corre- 
sponds to Rrim ~ 0.39 x (d/260pc) AU, or exactly half 
of the radius used in Model 1. In principle, this smaller 
radius will result in a dust temperature in the inner rim 
of ~2000K, using the opacity stipulated by the extinc- 
tion events. While the interferometry does suffer from 
a poor sampling of the u-v plane, this is a discrepancy 
which deserves a closer look. There are several ways of 



resolving the problem: the presented model includes a 
large cleared-out inner region due to the high tempera- 
tures of small grains. However, the location of the rim, 
given a dust sublimation temperature, actually depends 
on the dust opacity. A non-grey opacity will tend to ab- 
sorb more energy in the optical than it can emit in the 
near-infrared, increasing the temperature of the grains 
in the optically thin region of the puffed-up inner rim. 
This will tend to erode the inner edge of the rim, as 
small grains are heated above their sublimation temper- 
ature. This is reflected in the inner radius of Model 1. 
Conversely, larger grains with a more grey opacity will 
significantly decrease the temperature at a given radius 
of the inner rim, and consequently i? r i m will decrease for 
the same dust sublimation temperature. Specifically, us- 
ing the grain opacity for the outer (^2 AU) disk rather 
than the grains responsible for the UX Orionis events, 
wil l place the inner r im close to the value determined 
bv lEisner etaD (|2004Tl . i.e. 0.3-0.4 AU for an optically 
thick inner rim. This scenario then requires that small 
grains only appear in the upper layers of the puffed-up 
rim, while larger grains can penetrate close to the star in 
a flatter structure. This extra copmonent of large grains 
within the present inner rim does not strongly affect the 
model SED of VV Ser. 

To explore the possibilities offered by interferometry, 
we have calculated the visibilities for our model as well as 
several alternative structu res of the inn e r rim and com- 
pared them to those from Eisner et al. (2004). In Fig. 
El the structures of three different types of inner rim are 
sketched (Model 1 refers to the standard model discussed 
in the previous sections). Fig. [7] shows the calculated 
SEDs for the three models, while Fig. |H| shows the best 
fit visibilities from the three models (allowing the posi- 
tion and inclination angles to vary). In the following, 
position angles are those of the disk plane measured east 
of north. It is clear that given the position angle derived 
from the Spitzer images presented in paper II, the NW 
baseline visibilites can be well fitted, but the NS and SW 
baselines produce visibilities that are severely underpre- 
dicted by Model 1, i.e. the visibilities are o verresolved in 
the di rection of the disk plane. As noted bv lEisner et alJ 
(2004), the single point observed with the SW baseline 
may be significantly more uncertain than the statistical 
error bar indicates. We therefore assign little weight to 
this point. However, the bad fit to the NS baseline ap- 
pears significant. Given the sparse uv coverage, several 
widely different scenarios may provide better fits. 

First, assuming that the position angle of 10-20° from 
paper II is correct, the easiest way to amend the model to 
fit the near-infrared visibilities is by making the radius 
of the inner rim significantly smaller, thereby increas- 
ing the visibility along the major axis of the system, 
i^-band visibilities of a model in which the inner rim 
radius has artificially been placed at 0.4 AU has been 
calculated (Model 2 in Fig. |HJ|. A significant problem 
with this model is that it does not fit the SED because 
the small grains i n the inner rim reac h temperatures of 
2000 K. Note that lEisner et alJ l)2004|) finds a tempera- 
ture of 1500 K at this radius due to an assumption of 
a grey opacity for the inner rim. Additionally, a good 
fit to all three baselines can only be found for a posi- 
tion angle of ~ 165 ° and an inclination of 40°, consistent 
with that found bv lEisner et alJ l)2004fi . but inconsistent 



with the orientation of the su rrounding; nebulosi ty as de- 
scribed in paper II. Recently. Ilsella et al.l l)2006l) found a 
position angle of 60-120° using the same in terferometric 
data, roughly consistent with the result of lEisner et al.l 
( 2004) . If the single uncertain point of the SW baseline is 
ignored, a good fit can be found for the ~ 70° inclination 
and ~ 10° position angle found in paper II. These two 
possibilities need to be tested by further interferometric 
observations. 

Another option to reconcile the interferometry with the 
Spitzer images is a model in which large grains in the disk 
mid-plane reach smaller radii in an optically thin region 
within the outer puffed- up inner rim (Model 3 in Fig. |rj|). 
This structure is, in fact, quite plausible since the tem- 
perature can be maintained below 1500 K, as discussed 
above. This is presumably a natural scenario in a disk in 
which the large grains have settled to the mid-plane. We 
therefore construct a model in which the structure of the 
"standard" model is maintained (with small grains in a 
puffed- up inner rim), but with an additional component 
between 0.27 and 0.8 AU consisting of the large grains of 
the outer disk and with a structure following Eq. (J2J. 
Essentially, this model represents a rough simulation of 
the effect of grain settling and density-dependent evap- 
oration on the structure of the inner rim. This disk has 
basically two rims, both with a maximum dust temper- 
ature of 1500 K. The central idea of this scenario is that 
the SED is not strongly affected, since both dust compo- 
nents have a temperature of ~ 1500 K. At the same time, 
the visibilities can be fitted by a disk with the same posi- 
tion angle and inclination as that found in paper II using 
the morphology of the nebulosity seen in the IRAC and 
MIPS images of VV Ser. By tweaking additional param- 
eters such as the density structure of the grey dust and 
the exact location of the inner rim, an excellent fit to 
all the observables can be obtained. In particular, we 
find that a better fit to the visibilities is achieved if the 
surface density within 0.8 AU is lowered by a factor of 5 
to create an inner large grain component with roughly 
unity optical depth in the radial direction. This creates 
an inner rim that is not perfectly sharp, but radiates at 
a range of radii. This seems consistent with the struc- 
ture predicted by Ilsella fc Nattal (2005) using a density 
dependent grain evaporation temperature. 

Another option to fit the interferometer data is to add 
an additional point source offset by a few mas to produce 
additional visibility in the NW and SW baselines. We 
do not attempt to model this scenario here, but simply 
note that a few extra visibility data points will be able 
to distinguish between an additional point source and a 
small inner rim. 

Finally, it is possible that hot gas within the inner 
rim may produce a significant co ntribution to the K - 
band continuum as suggested bv lAkeson et al.l |2005), 
who showed models of PTI visibilities for a small sample 
of T Tauri stars. In fact, gas may produce much the same 
near-infrared morphology as that modeled here using two 
dust components. However, in order for the gas to be 
sufficiently luminous, a hi gh accretion rate of a t least 
lO- 7 M yr- 1 is required l|Muzerolle et al.ll2004jl . The 
accretion rate of VV Ser is lower than this, as evidenced 
by the lack of UV excess and the very low disk mass. 

6. CONCLUSIONS 
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Fig. 6. — Sketches showing the inner rim structure for three 
different disk models that have been compared to the interferom- 
etry and SED. Top panel: The "standard" model with the inner 
rim consisting of small grains at 1500 K, placing it at 0.8 AU. This 
model fits the SED, but not the visibilities. Middle p anel: Model 
using the inner rim radius of ~ 0.4 AU suggested bv lEisner et all 
(2004). Using the inclination and position angle from paper II as 
well as small grains in the rim, this model fits neither SED nor vis- 
ibilities. The visibilities can be fit better by changing the position 
angle and inclination significantly. Bottom panel: model using 
an inner rim at 1500 K for both grains populations present in the 
disk (small and grey grains). This model fits both the SED and 
visibilities with the position angle and inclination from paper II. 




0.1 1.0 10.0 100.0 

Wavelength [micron] 



Fig. 7. — Spectral energy distributions for VV Ser using three 
different models for the inner rim (see Fig. ©and Table|3J. 

We have presented Spitzer mid-infrared spectroscopy 
of the UX Orionis star VV Ser combined with additional 
spectroscopy and imaging data from the literature span- 
ning wavelengths from 0.1 to 850 /im. Using an axisym- 
metric dust radiative transfer model, we have reached 
the following conclusions: 
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Fig. 8. — Normalized _R"-band visibility amplitudes of the of the "standard" model (Model 1), the model with an inner r im at ,4 AU 
(Mode l 2) and the model simulating settling of larger grains to the disk midplane (Model 3). The observed visibilities are from lEisner et all 
(2004). The positition angle has been varied to find the best fit within -30 and 30°. The range was chosen to be reasonably consistent with 
the orientation of the bipolar nebulosity around VV Ser presented in paper II. The arrow indicates the position angle of the disk plane 
suggested by the Spitzer images presented in paper II (13°). The symbols indicate the observed uv points for the three baselines, while the 
curves show the model tracks (left panel) and visibilities as a function of hour angle (right panel). 



• All available data are consistent with the interpre- 
tation of the UX Orionis phenomenon in which 
the central star is transiently attenuated by dust 
clumps in the inner puffed-up region of an inclined, 
self-shadowed disk (the specific inclination is sup- 
ported by the imaging data presented in Paper 
II). While VV Ser is a very specific case that has 
been singled out due to the availability of a wealth 
of multi- wavelength, high quality data, it provides 
strong support for this interpretation in the more 
general case. 

• The disk appears to be both self-shadowed and low- 
mass, and both properties contribute to the quick 
decline with wavelength of the SED at mid- and 
far-infrared wavelengths. 

• The duration of the extinction events is consistent 



with an orbit of the dust perturbations in the inner 
rim at a radius of ~ 1 AU. 

• The dust responsible for the extinction events must 
be dominated by small, interstellar-like grains, 
while the dust in the outer (> 2 AU) disk ap- 
pears to be dominated by somewhat larger grains 
(~ 1 jum) . It is unlikely that the extincting dust 
clumps can be located significantly closer to the 
central star, since this would require them to sub- 
tend a much larger azimuthal angle than that cor- 
responding to the scale height of the inner rim. Ad- 
ditionally, they would be heated above their subli- 
mation temperature. 

• Our best-fit inner rim radius of 0.8 AU is in 
contradiction with the i-C-band interferometry of 
lEisner et alJ l)2004() which clearly indicates an in- 
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ner rim radius of 0.3-0.4 AU. A possible solution is 
to let the puffed- up inner rim consist of small grains 
at 0.8 AU while allowing a population consisting of 
larger grains to penetrate to 0.25 AU in a mid-plane 
that is optically thin to infrared photons. 

• The detailed structure of the inner rim described 
above is constrained both from the interferometry 
and the independent measurement of the disk incli- 
nation and position angles from large scale imaging 
as discussed in Paper II. 

• We the refore interpret the if -band interferometric 
data of lEisner et alJ (|2004) as evidence for larger 
grains in the disk mid-plane penetrating closer to 
the star than the small grains in the puffed-up inner 
rim. This is consistent with a scenario in which 
all dust grains sublimate at ~ 1500 K and larger 
grains have settled to the disk mid-plane, leaving 
only small grains in the disk surface. 

• Clearly, further tests of the detailed structure of the 
puffed-up inner rim are within reach of the current 
observational capabilities. Aperture synthesis ob- 
servations with upcoming facilities such as VLTI- 
AMBER in combination with multi-wavelength 
imaging and spectroscopy are certain to greatly en- 



hance our understanding of the innermost parts of 
proto-planetary disks. 
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